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ABSTRACT 

The three dimensional structure of the brightest part of the Orion Nebula is assessed in the light of 
published and newly established data. We find that the widely accepted model of a concave blister of 
ionized material needs to be altered in the southwest direction from the Trapezium, where we find that 
the Orion-S feature is a separate cloud of very optically thick molecules within the body of ionized 
gas, which is probably the location of the multiple embedded sources that produce the optical and 
molecular outflows that define the Orion-S star formation region. Evidence for this cloud comes from 
the presence of H2CO lines in absorption in the radio continuum and discrepancies in the extinction 
derived from radio-optical and optical only emission. We present an equilibrium Cloudy model of 
the Orion-S cloud, which successfully reproduces many observed properties of this feature, including 
the presence of gas-phase H2CO in absorption. We also report the discovery of an open-sided shell 
of [O III] surrounding the Trapezium stars, revealed through emission line ratio images and the onset 
of radiation shadows beyond some proplyds. We show that the observed properties of the shell are 
consistent with it being a stationary structure, produced by shock interactions between the ambient 
nebular gas and the high-velocity wind from 9^ Ori C. We examine the implications of the recently 
published evidence for a large blueshifted velocity of 9^ Ori C with respect to the Orion Molecular 
Cloud, which could mean that this star has only recently begun to photoionize the Orion Nebula. 
We show that current observations of the Nebula do not rule out such a possibility, so long as the 
ionization front has propagated into a pre-existing low-density region. In addition, a young age for 
the Nebula would help explain the presence of nearby proplyds with a short mass-loss timescale to 
photoablation. 

Subject headings: H II regions — ISM: individual: Orion Nebula, NGC1976 — star formation 



1. INTRODUCTION 

The brightest portion of NGC 1976, the Orion Nebula, 
is commonly called the Huygens region (after Christiaan 
Huygens, who published the first drawing of the nebula in 
1659) and its form in three dimensi ons was the su bject 
of many early papers (reviewed in O'Dell 2001a The 
presently accepted basic model of a photoiomzed thin 
layer of gas fiowing off the side of the Orion Molecular 
Cloud (OMC) facing us was invoked to explain the pro- 
gressiv e blueshift of the emission lines with r espect to the 
OMC ( [Zuckerman Il973[ [Balick et al.||1974[ ). The Bright 
Bar along the southeast border ot the Huygens region is 
a location where the main ionization front (MIF) is tilted 
almost along the line of sight to the observer. The other 
portions of the MIF have been mapp ed in three dimen- 
sions (3-D) by |Wen fc O'DeU] ( |1995| by a method that 
assumes all the ionizing radiation arises from the 07Vp 
star 9^ Ori C, and then calculates the distance to the MIF 
that would satisfy the known conditions of gas density 
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and observed Ha surface brightness, this being a broader 
applic ation of a method first presented in Baldwin et al. 
(1991 1. Although this 3-D map has the limitation that 
it becomes progressively less accurate as one moves away 
from the line of sight towards 9^ Ori C, it did establish 
that 9^ Ori C lies about^ 0.2 pc in front of the MIF, it 
confirmed the structure in the Bright Bar region, and 
showed that the nebula was otherwise a concave surface 
with a large bump to the southwest of 6*^ Ori C in the 
Orion-S star formation region. The surface brightness 
near 9^ Ori C can be explained as emission by a con- 
stant de nsity layer of about 0.1 pc thickness (Pogge et 
al.|1992 1 . This is only a reference number since the emis- 
sivity must be much higher near the MIF and drop as 
the square of the density. In a nebula where there is free 
expansion of the photoionized gas and there is a single 
dominant ionizing star, a concave shape of the MIF is the 
natural result, with bumps, valleys, and ridges refiecting 
underlying conditions of the host molecular cloud. The 
ionization front is expected to be closer to the ionizing 
star where the molecular cloud density is higher and fur- 
ther away where the underlying density is lower. The 
ionized gas density drops rapidly away from the MIF, 
but it is not clear what the density is in the immediate 
vicinity oi 9^ Ori C since it is expected that its intense 
high velocity wind would create a hot low density cavity 
around it. This cavity has not been detected directly ex- 

^ Throughout this paper we will adopt a distance of 440 pc to 
the Orion Nebula, a value derived | |0'DeIl fc Henney||2008| using 
the results of what are currently thought to be the best mdepen- 
dent determinations. Earlier papers that used different assumed 
distances will have their results scaled to this distance. 
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TABLE 1 

Heliocentric Velocities of the Ma,ior Orion Nebula Components 



Designation 

VoMC 

Vmif 

Vmil 

Vrrap 
^Cluster 
^fllOri C 
VArc 

^Veil-H+ 
Va l_A 

Vn i_B 
VH2C0 



Components VQ{kms ^) 

Molecules 25.8 ±1.7 

[Oi],[Sii] 25.5 ±1 
[Oii],[Oiii],[Nii],[Cliii],H+,He+ 18.0±1.4 

Ori A.e^ Ori B,el Ori D 24 ± 3 

Lower Mass Stars 25.8±1.0 

6li Ori C 13 

Big Arc 9.2 ± 1 

S++,P++,[Nii],Hei,[Oii] 3±2 

Hi 23.4 ±0.01 

Hi 19.3 ±0.03 

H2CO Absorption 28.6 ±1.6 



Source 



O'Del l et al. 

(2UUlb" 
(2001b 



20081 



'O Dell 
O'Dell 
Abt et al. ( 
^ uresz et al 
Stahl et aljlj^^ 
This study. 



Abel et al. 



Abel et al. 



(2006 
( 2DU5 



Abel et al. (20W| 

Manguni et al. 1n993|) 



Note. — Local Standard of Rest velocities can be obtained by subtracting 
18.1 km s^^ from Vq values. 



cept for observations of stan d-off shocks in front of the 
proplyds closest to 9^ Ori C ( iBally et al.||l998[ |2000 1. 

In addition to the rich Orion IN ebula Cluster (ONC) 
centered on the bright Trapezium stars, there are two 
star formation centers imbedded in the OMC. The first 



is associated with the deeply imbedded (0.2 pc, Doi et al. 

f04 ) BN-KL infrared and radio sources to the northwest 
0^ Ori C and the second in the Orion-S region. 
With the discovery of 21 cm absorption lines in the 
radio continuum spectrum (van der Werf & Goss 19891 
it was recognized that there was a Veil oi neutral ma- 
terial on the observer's side of the nebula and subse- 



2006 1 established its physical characteristics and approx- 



imate location of about one parsec on the observer's side 
(henceforth foreground) oi 9^ Ori C. Recent detaile d re- 
views have covered the ONC dMuench et al.||200 8D and 
the ne bula plus obscured star-formation regions ( O'D ell 
et al.lfeOOS, and see also, O'Dell 2001a, O'DcU 200TEl) 



The Huygens region occupies the northeast corner of 
a much lar ger structure called the extended Orion Neb- 



ula (EON, Giidel et al. 20081. It is known that there 



is a systemati c how of ionized ma terial into the EON 



(O'Dell 2001a 



20051, but the lower sur- 



Henney et al. 

t'ace brightnes s "Has limited the number o f investigations 
of this region ( Subrahmanyan et al.|200T ). However, the 
EON is the location of t wo X-ray bright regions of hot 
plasma (Giidel ct al. 2008) . High res olution opt ical ( Hci? 
[ney et al. 2007) and infrared (Megeath fc Robberto_2006j ) 
images of the EON are now available. Although these im- 
ages show many interesting large-scale features, there is 
neither the detail nor abundance of stars of the Huygens 
region. 

In this paper we will present the most relevant infor- 
mation about the ONC and the Orion Nebula, then inte- 
grate this into a modified picture of the Orion Nebula's 
3-D structure and its history. In § 2 we present the most 
useful information, then give the results of where this 
information leads in § 3. 

2. BACKGROUND INFORMATION 

2.1. Emission Line Images of the Orion Nebula 

The Huygens region has been the subject of numerous 
imaging studies. Arguably the most u seful groundbased 
study is that of |Pogge et al. (1992), which utilized a 
Fabry-Perot system to isolate emission from the Ha, H/3, 



quent absorption line spectroscopy (Abel et al. 2004b, l&'^oi 



[O III] 5007 A, [N ii] 6583 A and 6 548 A, [Sii] 6716 A and 
6731 A, and [Hei] 6678 A lines. lO'Dell fc Wong|(|l996 l 
presented a mosaic of Hubble Space T'elescope (HST'j 
WFPC2 images at the superior resolution of the HST. 
The particular advantage of the HST images is that the 
WFPC2 images allow clear discrimination of not only the 
isolated [O iii] 5007 A line, but also both the Ha 6563 A 
line and the nearby [N ii] 6583 A line in addition to the 
filters being narrow enough that they provide a good iso- 
lation of the emission lines against the strong scattered 
light continuum that primarily arises from dust grains in 
the dense region immediately beyo nd the MIF (O'DeU] 



1999p . A later HST survey ^Henney et al.||2007| ) 
with the ACS camera covered a wider field of view with 
pixels one half the angular size of those in the WFPC2, 
but, the filters used don't allow a clear delineation of the 



important emission lines (ODell 20041. The Huygens 
region has been mapped with the VEA at about 1.7" 
resolution ( [O'Dell fc Yusef-Zadeh]|2000[ ) at the extinc- 
tion free 20.5 cm conti nuum. By com paring the surface 
brightness calibrated (O'Dell & Doi' 1999) images and 
the radio images, it was possible to derive a map of the 
optical extinction across the Huygens region and to gen- 
erate extinction-corre cted versions of the HST WF PC2 
emission line images (O'Dell & Yusef-Zadeh 20001. In 
this study we employ the several forms of the WFPC2 
images, all processed in part using the IRAF package^. 
We see in Figure 1 that the Orion-S region to the south- 
west of the Trapezium resembles the Bright Bar region 
in being enhanced in low ionization [N ii] emission and 
being much brighter than adjacent areas. These char- 
acteristics are consistent with the brightest part of the 
central Huygens region being both closer to 9^ Ori C and 
also being an inclined face of ionized gas. This pict ure is 
consistent with the results of the 3-D modeling of |Wen| 
& 0'Dell|(p95|. 



2.2. The Orion Nebula's Veil 

Although the optical extinction of the Veil has been 
recognize d for some time, it has been est ablished only 



recently (O'Dell et al. 



1992 



O'Dell 2002 1 that most of 



this extinction arises from the near side of the ionized 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with the 
National Science foundation. 
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22.0 20.0 18.0 16.0 14.0 12.0 10.0 



Fig. 1.— This 233" x 219" image is composed of a mosaic of WFPC2 images ( |0'Dell fc Wong|1996| l with F502N ([O III]) as blue, F656N 
(Ha) as green, and F658N ([Nil]) as red. North is up and the labels along the edge depict the Higfit Ascension beyond 5:30:00 and the 
Declination south of —5:20:00 (2000). Major outflow features are labeled in addition to objects discussed in the text. Throughout this 
paper a position based designation is used (O'Dell &: Wen 1994} except for large individual features such as Herbig Haro objects. The 
white dashed ellipse encloses the smaller features collectively discussed in the text as the Orion-S feature. Red lines and circles represent 
strong infrared H2 features (Kaifu ct al. 2000 Stankc ct al. 2002) not associated with the BN-KL deeply embedded sources. The strong 
dark red/blue contoured lines indicate <JO outflow ^/japata et al.^,2005^ and the pas tel colored pink /light - blue contoured lines indicate 
SiO outflow ( Zapata et al.|2006 l. The orange contours are H2CO absorption features I Mangum et al.|l993 1 at the level of — 80K km s~^ 
(heavier) and — 50K km s~^ (lighter). The green contours of increasing thickness depict regions of the difference of extinction {c^jf as 
defined in § 2.4) of 0.1, 0.2, and 0.3. The sharp boundary of the more southerly c^jf excess feature is the result of reaching the edge of 
the WFPC2 field of view. The region most likely to contain the sources of the high velocity optical outflows is shown as a dark dashed 
ellipse (O'DoU & Honnoy 2008). Open squares indicate the positions of H2O maser sources (Gaume et al. 1998). Point sources within 
the dashed outline are coded by the shortest wavelength of their detection, with filled white squares indicating the positions of radio only 
vis ible sources (Zapata ct al. 2004a b, 2005), red squares the positions of sources seen only in the mid infrared (Smith et al. 2004 Robberto] 
|et al.|2005j , and ttlled orange circles positions of stars in t he near infrared catalog of .Hillenbrand fc Carpenter, ( ,2000 ). I'hc white contours 
show the 350 /xm emission in this area (Houde et al.|2004| in units of Jy per 12" bea!m. I'he irregular dashed white line indicates the field 
where Cjj/j was determined both by the radio/ optical and optical line ratio methods. 



4 



O'DELL ET AL. 



zone and not within it. The detailed analysis of Abel 



et al. (2005 20061 established that the physical condi- 
tions in the two 21 cm H i absorption velocity compo- 
nents of the Veil are rather different, the energy density 
in their Component A being dominated by the magnetic 
field measured from the Zeeman effect. The elements C, 
S, Mg, and Si are ionized in the Veil by the FUV (5- 
13.6 eV) radiation that penetrates it, even though the 
Veil is optically thick to EUV radiation. There must be 
a secondary hydrogen ionization front associated with 
the Veil and on the far side (a way from the obse rver and 
closer to 0^ Ori C) of the Veil. [Abel et al | ( |2006[ ) identify 
emission lines probably arising from the Veil's ionization 
front. 

The radial velocities of the various emission and ab- 
sorption line components along the line of sight through 
the Trapezium are summarized in Table 1. 

2.3. H2CO is Seen in Absorption in Orion-S 
An early study (Johnston et al. 19831 at 16" resolu- 



tion detected H2CO in absorption against the radio con- 
tinuum with the location identified as Orion-S , a region 



wi th mu ltiple known molecular emission lines (O'Dell et 
al. | 2008| and well define d bipolar r nolecular outfiows in 
CU (I Zapata et al.||2005'| and SiO ( |Zapat^^t_ar]|2006|) 
A hig her resolution (5.1 x 7.6") study ( jlVlangum et al. 
1993 1 confirmed the presence of this absorption feature 
and the results are shown in Figure 1 and Table 1. In 
Figure 1 we see that the H2CO absorption is distinct 
from the bipolar outflows, the strongest infrared and ra- 
dio sourc es, and the sources of t h e high velocity optical 
outflows ( [O'Dell fc Henney|2008] ). IMangum et~^ ( |1993 1 
point out that the presence of Ji2CO in absorption means 
that the cloud containing it must lie in front of the ion- 
ized gas, a posi t ion al so taken in more general form by 



Johnston et al.| ( |1983[ | and Wilson ct al. (2001). Since 
H2CO can only exist m a cold dense gas that is optically 
thick to FUV (and therefore also EUV) radiation, there 
will be a corresponding high extinction at visual wave- 
lengths. Therefore, we have looked for optical extinction 
associated with this feature with the results described in 
the next section. 

2.4. Anomalies in the Extinction in the Orion-S Region 

The optical appearance of the Huygens region is 
strongly affected by extinction occurring within the Veil. 
The clearest example is the Dark Bay feature to the east- 
northeast of the Trapezium. This extinction generally 
decreases away from the Dark Bay as the line of sight 
is moved to the southwest. The extinction has been de- 
rived in several fashions. In slit spectroscopy, the com- 
mon approach has been to compare the flux ratios of the 
strongest Balmer series lines with the ratios expected 
from theoretical predictions calculated for the local con- 
ditions (primarily the electro n temperature) . A go od ex- 
ample of this is the study of Baldwin et al. (1991), who 
obtained sample spectra along a well-defined east-west 
path beginning about 30" west of 9^ Ori C and derived 
the extinction from comparing Paschen lines with H7. 
The high est spatial resolution study is that of lO'Dell et' 
aI]( |2Ci03| , where cahbrated HST WFPC2 Ha and H/3 im- 
ages from a single WFPC2 pointing to the southwest of 
the Trapezium were compared with theory. Pixel by pixel 
extinction corrections were obtained as part of a study 



of electron temperature fluctuations. The widest field of 
view hi gh resolution determination o f the extinction was 
that of JO'Dell fc Yusef-Zadeh|[2000| ) who compared the 
Ha surface brightness of a mosaic of gaussian blurred 
HST WFPC2 images with VLA 20 cm images obtained 
at 1.7" resolution. In all these studies similar extinction 



curves were used (Costero fc Peimbert||1970 



Cardelh et 



al.|1989 1 . The derived extinction is conimonly expressed 



as ch/3, the logarithm of the extinction at 4861 A, the 
wavelength of H/3. In this work we will refer to the ch/3 
value derived from the comparison of the radio contin- 
uum and Ha as c^^^and to the value derived from com- 
parison of WFPC2 Balmer Ha and H/3 images as c"°, 
while the difference of the two will be c^^^ = c^^^ — c^J^ . 

Although a comparison of the results for c^J^* and c^J^ 
for nine well-sampled smaller regions showed that there 
was a good agreement of the results o f the two meth- 
ods for finding ch/3 (ratio of 0.93±0.19,( |O'Dell fc Yuse? 
Zadeh|200 0|), there is reason to expect that the two meth- 



ods should not agree exactly. 

The first reason would be the effect of scattering of the 
optical Hues by dust in the dense photon dominated re- 
gion (PDR) that lies immediately beyond the MIF. This 
scattering is what produces a broad redshifted compo- 
nent in the high resolution spectra of the intrinsically 
narrower heavy ion emission lines and it must also be 
present in the Balmer lines. This is why the continuum 
of the Orion Nebula is mu ch stronger than exp ected from 
simply atomic processes (Baldwin et al. 1991), in effect. 



the Orion Nebula is also a bright reflection nebula. One 
would expect the effect of this process to cause an un- 



derestimate of c^J^* 



and would affect c^^ less since one 



is dealing with only the difference in the amount of scat- 
tering of the Ha and H/3 lines, rather than the absolute 
amount at Ha, as done in deriving c^J^*. Taking the frac- 
tion of the light scattered in the forbidden lines, this po- 
sition dependent effect should only be a few hundredths 
in the logarithm (dex). 

The second reason to expect c^Jf to be non-zero is that 
when using only the optical emission (the Balmer line 
ratio method) one is dealing entirely with radiation that 
has not suffered a large amount of extinction, whereas the 
radio continuum to Ha method can be thrown off by the 
fact that some of the total radio continuum is produced 
by volumes behind regions that are very optically thick 
to Ha. In this case, one would be using a radio signal 
from all along the line of sight but Ha emission from only 
in front of the obscuring material and c^J;;* would be too 
large, i.e. c°Jf would be positive. 

The third reason for a difference would be variations 
in the ratio of total to selective extinction. Derivation 
of depends upon the reddening curve adopted, which 
is primarily determined by the ratio of total to selective 
extinction, while c^J^* is less dependent on this ratio. 

Since the presence of H2CO in absorption indicates 
that near Orion-S there is a region where radio contin- 
uum emission occurs in part behind the totally obscuring 
cloud, we looked for signiflcant variations in c°Jf in the 
region southwest of the Trapezium where this could be 
done. The results are also shown in Figure 1, where we 
have superimposed the contours of large values of c^]^. 
One sees that there are two regions where c™f is large, 
the northerly one nearly coinciding with the H2CO ab- 
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sorption and the other coinciding with the distinct ob- 
scurinR feature to its south d esignated as the SW Cloud 
( [O'Dell fc Yusef-Zadeh|pOOO| . This means that in both 
these regions optical emission lines are recording only 
part of the ionized gas along the lines of sight. Another 
way of wording this is that in these regions the model 
of the MIF front based on a sing le dominant emitting 
layer cannot be accurate. Since the 'Wen & O'Dell ( 1995 1 
study used the reddening corrected Ha line for reference, 
in these regions the model will be a depiction of the closer 
(to the observer) emitting region. This means that the 
"bump" in the Orion-S region is there, but it is not sim- 
ply the result of a nearer part of the MIF, but, it is an 
isolated optically thick feature directly ionized on its sur- 
face by 0^ Ori C and behind it is another ionized region. 
The more southerly region of large c^^^^ is probably sim- 
ilar, but it must not be as dense as the Orion-S feature 
since one does not see H2CO in absorption. 

2.5. Radiation Shadows in the Orion Nebula 

Nature provides a way of probing where material is 
located within the Orion Nebula through the ionization 
shadows created beyond the proplyds. In a study of lin- 
ear features withi n the Orion Nebula and the Helix Neb- 
ula, [O'DcU (2OOOI) established that these features can oc- 
cur when there is a single dominant ionizing star and 
an intervening object optically thick to the EUV Lyman 
continuum (LyC) radiation. The first orde r theory for 
this s ituation had already been determined ( Canto et al. 
1998| . Gaseous material lying beyond an object like a 
proplyd (in the Orion Nebula) or a knot (in the Helix 
Nebula) will be shadowed from photoionizing LyC pho- 
tons coming directly from the ionizing star, but will be 
exposed to LyC photons produced through the recom- 
bination of surrounding nebular gas. This diffuse radia- 
tion field will be a small fraction of the direct radiation 
from the ionizing star and there will be fewer higher en- 
ergy photons than in the direct radiation. The differ- 
ence in the radiation field is because the recombinations 
occur preferentially to energy states just above hydro- 
gen's ionization threshold. The lower density of ionizing 
photons means a lower rate of photoionization and the 
different energy distribution means that any shadowed 
gas will equilibrate at a lower electron temperature. The 
lower temperature means that forbidden line radiation 
will have a reduced emissivity and recombination lines 
will have an increased emissivity. This shows up partic- 
ularly well when looking at the ratio of i mages in [O iii 



(20001 



and Ha, as shown in Figure 2 and Figure 3. O'Dell 

established that the linear rays seen in the [O mJJ/- 
age occurred along lines of sight passing through 6^ Ori C 
and a proplyd, but not all proplyds have visible radiation 
shadows. 

These shadowed regions will only appear when there 
is material present. This means that the absence of a 
radiation shadow indicates that only low density gas lies 
within that conic column. In the case where the radia- 
tion shadow does not begin at the proplyd, then the gas 
occupying the space between the proplyd and the start 
of the radiation shadow must be of low density. 

Figure 2 has had all the linear and nearly linear fea- 
tures traced, solid lines being used for the indistinguish- 
ably linear features and dashed lines for the nearly linear 
features (which cannot be radiation shadows unless sec- 




FlG. 2. — The same field of view as Figure 1 is shown exeept 
that here the ratio of the F502N over F656N images are shown. 
Solid lines depict radial ionization shadow features that project 
back towards 9^ Ori C and dotted lines depict significant nearly 
linear features that do not project back to 6^ Ori C. Almost all 
of the solid-line features have an identified proplyd which is the 
occulting source and in this case the designation of the proplyd is 
given, adjacent to the proplyd's position. The exception to this 
rule is the ionization shadow feature lying beyond the HH 202-S 
shock. The white dashed line shows the boundary where radial 
shadows begin. 

ond order effects discussed in 'O'DcU 2000 are in play) . 
One sees that in essentially every case of a linear fea- 
ture there is an identified shadowing proplyd and, again, 
that not all proplyds have radiation shadows. In many 
cases, the associated proplyd is well inside the radia- 
tion shadow (closer to 0^ Ori C). The general pattern 
is clear, none of the proplyds near 0^ Ori C have radia- 
tion shadows extending close to them and the proplyds 
with closely approaching radiation shadows are distant 
from 6^ Ori C. We interpret this to mean that there is 
a region of low density (or very different ambient condi- 
tion) gas lying within the approximate boundary shown 
in Figure 2. This irregular boundary in the plane of the 
sky is the projection along the line of sight and the cen- 
tral cavity is probably somewhat larger than the outlined 
boundary. The irregular boundary drawn indicates that 
there is an asymmetric distribution about the dominant 
ionizing star 0^ Ori C. 
It sh ould be noted that the h igh velocity feature HH 

also casts a radiation 



202-S d O'Dell fc Hemiey 
'i'his indicates t 



20081 

lat this shock produces a suf- 



shadow. 

ficiently high density that it is optically thick to LyC 
photons and becomes a local ionization front within a 
generally ionized volume. This vitiates any argument 
that HH 202 was caused by high velocity material from 
the Orion-S region impinging on cither the foreground 
Veil or a section of the MIF that had curved towards the 
observer and intercepted the coUimated outflow. 

2.6. Large Scale Ionization Structures within the inner 
Huygens Region 

Within the paradigm of the Orion Nebula being de- 
scribed to first order as a relatively thin blister of ion- 



6 



O'DELL ET AL. 



ized gas on the 6^ Ori C side of the MIF, one expects 
that the [N ii] emitting zone will be quite narrow as com- 
pared with the Hof emitting zone, which includes all of 
the ionized gas. [O iii] emission must arise from closer 
to 9^ Ori C and over a larger volume since the effective 
temperature of 9^ Ori C is insufficient to produce a He++ 
zone (where 0++ has been ionized and [O iii] is not emit- 
ted). This means that the [Nil] /Ha ratio is a good in- 
dicator of when the MIF is highly tilted with respect to 
the plane of the sky. This is the reason why the Bright 
Bar feature is so well delineated in the left panel of Fig- 
ure 3. Although there are many unexplained features of 
the Bright Bar, the overwhelming evidence is that it is 
an escarpment in the MIF where we see almost along the 
local ionization front. This is probably also the correct 
interpretation of the low ionization zone to the south- 
west of the Trapezium that lies on the northeast side of 
Orion-S. 

2.6.1. A Previously Unrecognized High Ionization Structure 
within the Inner Huygens Region 

Figure 3 (left panel) reveals an unexpected large scale 
structure in a C-shaped zone of low [Nil] /Ha with 
9^ Ori C located near its middle. Close examination of 
Figure 3 (right panel) indicates that where [Nil] emis- 
sion is weak, [O iii] emission is strong. This means that 
the C shaped zone is an ionization phenomenon. The 
fact that the [N ii] deficit and [O iii] excess is an irreg- 
ular line, rather than a closed area, indicates that this 
is a shell of high ionization. Henceforth we will call this 
the [O III] Shell, with full recognition that this feature is 
probably not closed and we do not understand its three 
dimensional characteristics. South of the Trapezium and 
Orion-S lies a large high ionization feature labeled as 
Big Arc South and Big Arc East. This was originally 
discovered in a low velocity resolution study (O'Dell et 



al. 1997), then included in a later higher resolution study 
" pJoi et al. 2004) , and the correct interpretation as a par- 
tial shell blueshifted about Vp) = 9 km s~^ being pre- 
sented in Garci'a-Di'az et al. ('2008') and confirmed by our 
remeasuring a sample ot the Bi g Arc from the [O iii] spec- 



tra data set (Doi et al 



2004). Although the Big Arc 



is seen in projection to be almost continuous with the 
southern arm of the [O iii] shell, it is very distinct kine- 
matically, and the physical relation between these two 
features remains unclear. 

The approximate alignment of the boundaries of the 
radiation shadows and the [O iii] Shell argues that this 
is an incomplete shell with lower central density within 
the larger mass of ionized gas. It is not an ionization 
boundary because one does not see [Nil] enhancement 
on its outer perimeter. The feature seems to open to 
the west-southwest, being most open in the direction of 
Orion-S. There is a nearly linear border to the high ion- 
ization zone, which we have labeled in Figure 3 as the 
SE-NW Ionization Boundary. It is not immediately obvi- 
ous that this feature is an ionization boundary or caused 
by the higher extinction of the foreground Orion-S Cloud 
(cf. § 3.1) 

Figure 3 also shows the boundary of the radiation 
shadows from Figure 2 and we see that there is good 
agreement of this boundary with the [O iii] Shell, which 
further strengthens the argument that this [O iii] excess 
traces the outer boundary of a volume of low gas density. 



We also show in Figure 3 the location of proplyds with 
stand-off shocks and find that all of them fall within the 
C-shaped [Oiii] Shell. Since it is most likely that these 
stand-off shocks are formed by direct exposure to the 
high velocity stellar wind oi 9^ Ori C, this argues that 
the [O III] Shell is created by that stellar wind, an idea 
developed in § 3.1. 

2.6.2. Physical Conditions in the [O in] shell 

The average hydrogen number density of the [O iii] 
shell, n, can be estimated from the Ha surface-brightness 
(photons cm^ s~^ sr~-^), given by 



^(Ha) = X cCt cfin^ L / ^t: , 



(1) 



where L is the path length of the line of sight through the 
shell and Xc is the electron fraction. Both hydrogen and 
helium are assumed to be singly ionized, with He/H = 
0.0977 ( |Esteban et al] p98), so that x^ = 1.0977. The 
effective line recombination coefficient aoff has a very 
weak density dependence and can be appr oximated in 



the Case B limit as 1.16x 10-^3T4"^ cm^ s"^ ([Osterbrock 



fc Ferland||2006| , where = T/IO'' K. For a spherical 
shell of radius R and fractional thickness A ^ 1, the 
maximum path length through the shell is L ~ 2A^/^i?. 
The observed radius and thickness of the [O iii] shell vary 
with PA by about 50% but we will take i? = 2 x 10^^ cm, 
A = 0.28 as typical of the well-defined northern arm of 
the shell, so that L ~ i?. 

The surface brightness of the shell is most easily es- 
timated from the sharp brightness jump at the inner 
edge of i ts limb-brightened nor thern arm. Using WFPC2 
images (lO'Dell fc Wong 19961, flux-calibrated following 
O'Dell fc Doi J 19991) and ext inction-corrected following 
O'bell k Vusef'-I<!!acieh| ([2000)), we find S'(Ha) = 2.85 x 



10^" photons cm'^ s ^ sr ^, which represents about 25% 
of the total nebular surface brightness at the shell posi- 
tion. We also find [Oiii]/Ha = 1.62 and no evidence of 
any [Nil] emission from the shell ([Nii]/Ha < 0.0 1). 

For the shel l temperature, we use the results of |Mesa- 
Delgado etal.| f2008 ) who find T ~ 8300 K at the position 
of the shell from the ratio of the 4363 A and 5007 A [O iii] 
lines (upper-right panel of their Fig. 6). This tempera- 
ture represents an average of the emission from the shell 
and the background nebula at that position, but there 
is no evidence for any significant jump in T at the inner 
edge of the shell, so we can be confident in using it for 
the shell temperature, giving an effective recombination 
coefficient of 1.36 x 10" ^■^cm^ s"^. 

Combining all the above, we find a shell emission mea- 
sure of EM — Xcn^L — 2.56 x 10^^ cm"^ and hence a 
density n — 3400 cm""^. The ionizing flux at the position 
of the shell is $h = Qb_/^t:R^ ~ 3.6 x lO^^ ^m-^ s-\ 
where Qh is the ionizing lu minosity of 9^ O ri C, assumed 
to be 1.8 X lO''^ s-i (Henney et al.||2005a |. 

The ionization parameter in the shell (<1>h/'^c = 0.35) 
is very large when compared with typical values for the 
nebula as a whole (0.01-0.02). Therefore, if dust grains 
are not significantly underabundant, the grains will dom- 
inate the ultraviolet opacity of the shell. Assuming a dust 
extinction cross-secti on of tJd = 10~^^ cm^ per H nucleon 
(Baldwin et al. 1991), the optical depth of grains in the 
shell to ionizing radiation is na^R^ — 0.2. For compari- 
son, the neutral hydrogen optical depth is much smaller: 
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Fig. 3. — The same field of view as in Figures 1 and 2 is shown in both panels. The left panel shows the ratio of F658N/F656N ([Nll]/Ha) 
and the right panel F502N/F656N ([Olll]/Ha). Both are intended to show the systematic ionization changes across the Huygens region. 
The left panel is little affected by local variations in interstellar reddening, while the right panel is. The dashed line indicates the boundary 
where radial shadows begin. The inner boundary of the [Olll] Shell feature discussed in § 2.6 is outlined with a curved solid line. The 
lighter solid line indicates the southwest boundary of a high ionization region close to the Trapezium stars. The circles are centered on the 
proplyds that have stand-off shocks, which indicate that they lie within the unshocked hypersonic stellar wind from 6^ Ori C. 



aBEM/$H ^ 0.022. Thus, the shell is optically thin to 
EUV radiation, which, combined with the high ioniza- 
tion parameter, means that low-ionization lines such as 
[N ii] are expected to be very weak. Indeed, no detectable 
[N ii] emission is observed from the shell. 

2.7. Space Motion of 6^ Ori C 

In spite of its brightness, it has been difficult to deter- 
mine the three-dimensional space motion of 9^ Ori C, the 
dominant ionizing star. This is true for both the tangen- 
tial and radial velocity components that are needed for 
deriving this motion. 

The brightness relative to the other cluster members 
presents particular challenges in determining the tangen- 
tial velocity. The most accurate study is that of vaii] 

(1988), who determined that 6*^ Ori C was 

-1 



Altena et al. 



moving at 4.8 ±0.5 km s towards Position Angle 142° 
and that this was significantly larger than the disper- 
sion value of 1.5 ± 0.1 km s^^ found for the ot her clus- 
ter star s (49 cluster members in the sample). [Jones& 



Walker ( 1988} did not include the brightest members of 
the ON C in their astrometric study, but did find a veloc- 
ity dispersion of 1.5 ± 0.7 km s""'^ for the brighter stars 
in their study. Tan (20041 interpreted this high tan- 



gential velocity away from the direction of the imbedded 
BN-KL cluster stars as evidence that 0^ Ori C shares a 
point of origin with the radio source BN, which is moving 
in the opposite direction at a high velocity (38 km s^^ 
at PA = 322°), and that the present tangential velocity 
would have placed the BN object in the vicinity of the 
Trapeziu m about 4 000 years a go. However, subsequent 
studies ( Rodriguez et al. 200"5|jG6me z et a l.pOOS, 2008] 
of the radio sources within the BN-KL region show that 
three radio sources in the center of the BN-KL region are 
moving away from a common center with a timescale of 



50 years and that the BN-KL source I used fo r reference 
by TmT (20041 was one of the moving objects. [Gomez et| 
jal^(2005) found that BN was moving towards the north- 
west at 24 km s^^ in the Orion cluster rest frame. The 
separation of 0^ Ori C and the dynamic center of BN- 
KL is 63", which means that a backwards p rojection of 
9^ Ori C's tangential motion as determined by van Altena 



et al. 



( |1988) would cross that position in about 28,000 
years. Even if it is moving away from the BN-KL center, 
this timescale indicates that it is unlikely 9^ Ori C is a 
runaway object from the BN-KL region event that pro- 
duced the three escaping radio sources and t he bipolar 
outflo w with an upper limit age of 1000 years (Doi et al. 



2002 1 that has produced the well known H2 fingers 
Even the modest anomalous tangent ial velocity of 4.8 



km s may be too high. The data of van Altena et al. 
(|1988) indicate that 9^ Ori C and 9^0m B are moving m 
opposite directions at about 3.5 mas per year. In a com- 
pilation of visual astrometric measur ements of the rela- 
tive positions of the Trapezium stars, Allen et al. ( 1974 1 
present results for these two stars over a time interval of 
125 years, during which time the stars would have sepa- 
rated by 0.43", but, there is no detectable change to an 
accuracy of about 0.15". 
Determination of the radial velocity of 6*^ Ori C is also 



diflflcult (O'Dell 2001b Vitrichenko 2002 1 . In this case 



the brightness of the star aids the collection of high ve- 
locity resolution spectra, but the star is complex, being 
an oblique-rotation magnetic star which shows a bewil- 
dering pattern of periodic and erratic spectral changes. 
It has a resolved early type star companion. Proper- 
itiesof the system are summarized in Mue nch et al.| 
20081). Utilizing on ly spectroscopic radial velocity data, 



itrichenko ( 2002 1 determined a systemic radial velocity 
of 10.9 ± 2 km s~ and by combining spectroscopic and 
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astrometric ( Kraus et al. 2007 Patience et al. 2008 ) data, 
Stahl et al. ( |2008 ) have determined the systemic radial 



velocity of 6 Ori C as 13 km s~ , which is quite differ 
ent from that of the other Trapezium stars and the lower 
mass ONC members, as summarized in Table 1, and indi- 
cating that it is moving rapidly away from the OMC and 
its host star cluster (relative velocity about 13 km s~^). 
In § |3.4| below, we examine the effect that such a rapid 
motion for 9^ Ori C would have on the present structure 
and evolution of the Orion Nebula. 

This special characteristi c (high spatia l velocity) may 
not be all that unexpected. Stone ( 1991| ) concludes that 
in general 40% of all O stars are runaways, although 
only one in ten OB stars has a known binary compan- 
ion. There is an observational selection effect operating 
in the ONC that favors detection of blueshifted radial ve- 
locities. A redshifted runaway would be moving into the 
obscured portion of the OMC and we would not be see- 
ing the photoionized Orion Nebula at visual wavelengths. 
This means that if we are going to detect an anomalous 
radial velocity, it is going to be a blueshift. The differ- 
ence in radial velocity with respect to the other Trapez- 
ium stars would mean that either this juxtaposition of 
massive stars is coincidental or perhaps that dynamical 
interactions with the other stars has produced 9^ Ori C's 
anomalous radial velocity. 

Certainly, runawa y stars are expected from a grouping 
like the Trapezium (Pflamm-Altenburg & Kroupa|2006l 
and have been created previously in the ONC. The diver- 
gence from Orion of the two high velocities stars AE Aur 
and /i Col has be en known for over a half century ( Blaauw 
fc Morgan 1954 ). Calculations of the t rajectories of these 
stars and t Ori ( Hoogerwerf et al. 12001 1 indicate that they 
all originated about 2.5 million years ago from a position 
now occupied by the ONC. This is not to argue that 
these stars arose in the Trapezium grouping, however, it 
is evidence that there was an earlier epoch of massive 
star formation in this vicinity. Given today's presence of 
three centers of star formation (the ONC, BN-KL, and 
Orion-S), it is not hard to accept that there was at least 
one earlier similar center. Likewise, the earlier creation 
of runaways argues that this could have again occurred 
and 9^ Ori C is the product. 

3. DISCUSSION 

The new features presented in § 2 have caused us to re- 
assess several major properties of the Orion Nebula. In § 
3.1 we consider the effects of the high velocity wind from 
9^ Ori C. In § 3.2 we present arguments for the Orion-S 
feature being a separate cloud within the ionized cavity 
of the nebula and in § 3.3 present a model consistent with 
its producing H2CO in absorption. We conclude with a 
discussion of how our view of the nebula would be altered 
if 6*^ Ori C has a large blueshifted velocity with respect 
to the other members of the ONC. 

3.1. The Effects of the High- Velocity Wind from 
9^ On C 

The presence of a high velocity stellar wind from 
9^ Ori C is characteristic of early spectral type stars and 
will certainly affect the conditions of the nearby nebular 
gas. In this section we address the nature of the wind, 
its expected interactions, and compare the expectations 
with what is observed. 



The radiatively driven stellar win d is modified by th e 

2002t, 



star's strong dipolar magnetic field ( Donati et al 
lea ding to a very anisotropic outflow ( [Smith ^ 



uller- 



|ton|| 2005| . However, the large inclination between the 
magnetic and rotation axes ( jWade et al. 2006 1 means 
that much of this anisotropy is smoothed out over the 
15.422 day rotational period ( [Stahl et al.|1993j|1996[ ), so 
that it is reasonable to assume that the wind is approxi- 
mately isotropic when considering its interaction with the 
Orion Nebula. Angle-averaged values for the mass-loss 
rate and terminal velocity of Af^ — 4x 10~^ Mq yr~^ and 
Kv = 1400 km s^^ have been determined from m odeling 
the X-ray emission from the base of the wind (Gagne 
eFar][2005| . 



A totally independent measurement of the strength 
of the fast stellar wind is provided by the [Oiii]/MIR 
arcs seen in front of the closest proplyds to 9^ Ori C 
playward et al.[|1994l iBally et al.[[1998[ [Robberto et al.| 
|2D0 5 Smith et al. 12005 ) , which are interpreted as the 
■ ■ otf 



stationary stand-on bowshocks that result from the in- 
teraction of t he proplyd photoevaporation flo w with the 
stellar wind ( Garcia- Arredondo et al. 20011. The po- 
sition of these bowshocks is determined by the balance 
between the ram pressures of the proplyd outflow and 
the fast wind from 9^ Ori C, so that they can be used 
as probes of MwKv so long as the proplyd parameters 
can be determined with sufficient pre cision. Application 
to the best- s tudied proplyd 167-3 17 (Garcia- Arredondo 



et al.[|200l[ [2002| [Henney et"al][2002 ) yields a value 



of (M^/10-^ Mq yr-^)(V^ /mO km s 



4.6 ± 1.4, 



which is consistent with the Gagne et al. (20051 values 
given above. 

3.1.1. Theoretical Predictions of the Wind/Nebula 
Interaction 

The interaction of the fast wind from 9^ Ori C with 
its surroundings differs from the standard spherical 
hot-shocked wind bubble scenario (see e.g., Dyson & 
Williams 1997) due to the effects of geometry, density 
gradients, mass- loading from embedded sources, and pos- 
sibly thermal conduction. The different wind interaction 
regions are summarized in Figure 4, which synthesizes 
findings from many recent theoretical and observational 
studies. Each region is discussed further below, starting 
from the inside. We concentrate on the inner parts of 
the nebular that are relevant to the observational mate- 
rial presented in this paper, leaving aside considerations 
of the EON from which x-ray evidenc e of stellar wind in - 
teraction has recently been reported ( Giidel et al.|2008 l. 

Mass loading zone — Mass loading of a stellar wind by 
embedded sources lowers the stellar wind velocity of the 
unshocked wind and consequently reduces the temper- 
ature and modifies the cooling properti es in the hot 
shocked wind region ( Hartquist et al.|1986 1 . In the Orion 
Nebula, the sources of mass are the photoevaporated 
fiows from the close-in proplyds (PB in the figure) and 
the stellar winds from the other Trapezium stars (SB). 
Of these sources, the proplyds are the most important 
( Garcia- Arredondo et al.[[2001 l, since the mass-los s rate 
from each proplyd i s of the same order as Mw (Hen 



ney fc 0'Dell|[1999[ [Henney et ar][2002| ), whereaslEe 
most massive of the other Trapezium stars are early B- 
type, with winds at least an order of magnitude weaker 
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Fig. 4. — Schematic cross-section through the inner Orion Nebula, showing the different zones of stellar wind interaction. Blue 
indicates stellar wind material, red indicates photoionized nebular material, and green indicates atomic/molecular gas, with darker shading 
corresponding to qualitatively higher density within each region. Arrows show the approximate direction of gas flows, with sizes that 
roughly scale with the momentum flux. Shocks in the stellar wind are shown by short-long-dashed lines, while the contact discontinuities 
are shown by dotted lines. The drawing does not represent any particular cut through the nebula, and is not exactly to scale, although the 
sizes in parsecs of some major features are marked. The observer is located off the bottom of the picture. The inset box shows a zoom of 
the central region. 

than that of 6*^ Ori C (MS in the figure). Unfortunately, 
the efficiency of the mixing between the stellar wind and 
the proplyd pho toevaporation flows is poorly constrained 
( Garcia- Arredondo et al._,2002| ), leading to large uncer- 
taities in the degree of deceleration of the wind. However, 
as long as the mixed wind remains highly supersonic, its 
momentum M^Vw is conserved. 

Shocked wind — When this mass-loaded, but still su- 
personic, wind interacts with the nebula, a two-shock 



flow p attern is formed ( Pikerner|1968 Dyson & de Vries 
1972|): an outer shock that accelerates the nebular gas 



and an inner shock (blue double-dashed line in Figure 
4) that decelerates the stellar wind. The two shocked 
regions are separated by a contact discontinuity (blue 
dotted line), across which the pressure is continuous. If 
the inner wind shock is non-radiative (as is usually the 
case), the shocked stellar wind (SW) will have a temper- 
ature T/IO^ K ~ 0.15(Kv/100 km s"^)^ and a thermal 
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pressure roughly equal to the ram pressure at the inner 
shock. 

Large-scale density gradients in the adjoining molecu- 
lar cloud (MC) mean that the inner shock is much closer 
to the star in the cloud direction (upper-left in Figure 
4) than in the opposing direction (ofF-picture to lower- 
right), leading to a higher pressure of the shocked stellar 
wind in the cloud-facing direction, which ther efore flows 
towards the low-density end of the bubble (Arthur & 
Hoarel'2006) in a transonic wind exhaust (WE in l* 'igure 



4). A simple steady-state model then gives the radius of 
the wind shock as 0.73 times the radius of the contact 
discontinuity. 

Shocked nebula — In the early evolution of an Hii re- 
gion, the ionization front is likely to be trapped inside 
the de nse shell swe pt up by the outer wind-driven shoc k 



Garcfa-Segura & Franco 1996 Arthur & Hoare 20061 



The density in the shell decreases with time, so that 
the ionization front eventally escapes from the shell and 
drives a low- vel ocity neutral shock ahead of itself into the 
molecular gas (IArthurP2007|. After this time, the outer 
wind-driven shock will be interacting with photoionized 
nebular g 3tS, clS shown in Figure 4. If this gas is static, 
with a constant hydrogen number density n cm~'^, then 
the outer shock will degener ate into a sound wave af - 
ter ~ 1.3/n^/^ million years (Dyson & Williams 1997), 
leading to dissipation of the shocked nebular shell, after 
which one will have approximate static pressure equilib- 
rium between the H ii region and the shocked wind. In 
the cloud-facing direction, where the ionized density is 
n = lO'^-lO* cm~'^, this will occur on a timescale much 
shorter than the age of the nebula, so that stagnant zones 
of ionized gas (SI in Figu re 4) can be confined by the wind 
( Arthur fc Hoare]|2006 1 . In the direction away from the 
cloud, the ionized density is much lower, n < 100 cm~'^ 
( Felli et ar]|1993 1, so that one may still have a dynamic 



shell (US) of swept-up H ii region that moves outward at 
a few tens of km s~ . 

Additional structures will arise due to the internal dy- 
namics of the ionized gas in the H ii region. Turbulence 



and gravita tional collapse in molecular clouds ( Vazquez- 
Semadeni et al. 2005) produce density concentrations on 



multiple scales (idealized as dense cores, DC, in Figure 
4), which retard the progress of the ionization front. In 
regions where the ionization front is locally convex or flat, 
photocvap oration flows (PF ) will develop (Henney et al. 
[2b05a; ^Mellema et al. 2006 1 , whi ch reach mildly super- 



sonic velocities of 20-40 km s (Dyson 



1968 



Bertoldi 



19891. Where one of these photocvaporation flows inter 



acts with the stellar wind bubble, a stationary shell (SS) 
of shocked ionized gas will form at a position where the 
ram pressure of the photocvaporation flow balances the 
thermal pressure of the shocked wind. 

3.1.2. Observational Evidence for the Ejfects of the Wind 
on the Ionized Nebula 

As reported in § 2.6.1 above, we have found an in- 
complete high-ionization shell in the inner Orion Neb- 
ula, which may be related to the action of the stellar 
wind from 0^ Ori C. Since the shell has both an inner 
and outer boundary (Figure 3), it would correspond to 
a stationary shell in the context of the above discussion 
(SS in Figure 4) . This means that the kinematics of the 



shell is likely to be unremarkable, as compared with the 
kinematics of the general [O iii] emission from the neb- 
ula. Just as with the rest of the nebular gas, the gas 
in the shell will flnd it easier to move in the direction 
away from the molecular cloud and so should be slightly 
blue-shifted on average compared with the systemic ve- 
locity of the stars and molecular gas. In spection of the 
[O III] velocity cubes of Doi et al. ( 2004 1 shows that the 
line proflles from the shell are, indeed, very typical of the 
nebula as a whole. 

The inner boundary of the shell corresponds to the con- 
tact discontinuity with the shocked wind bubble, so the 
shell thermal pressure should equal that of the shocked 
wind. From the p hysical parameters derived from obser- 
vations in 
(1 -I- Xc)nT 



2.6.2 



one finds a shell pressure of Psh/k 
X 10^ cm~^ K. The bubble pressure 
will be equal to the wind ram pressure at the inner wind 
shock: Pbub/fc = MwKv/47ri?2^fc ~ 7.84 x 10^ cm'^ K, 
where we have assumed that the inner shock radius is 
i?in — 0.73i? (see above) and that the wind momentum 
loss rate M^Kv is the proplyd-derived value discussed 
in § 3.1. Given the uncertainties in their derivation, 
these pressures are indistinguishably the same, providing 
strong evidence that the [O iii] shell has been compressed 
by the stellar wi nd. 

The Red Bay ( IGarcfa-Dfaz fc Henney|2007| IS a region 
in which the high-ionization g traced by [Siii] and 
[Oiii] lines, has a velocity similar to that of the molec- 
ular cloud {Vq — 25-30 km s^ ), rather than the more 
blueshifted velocities {Vq < 25 km s^^) found elsewhere 
in the nebula. The region of relatively redshifted veloci- 
ties starts at the Trapezium and extends in a broad strip 
(~ 75" X 25" ) to the east-southeast, which can be best 
appreciated as a light-co lored region in the fO iii ] mean 



velocity map (Fig. 13 of Garci'a-Di'az et al 



2008|.« The 

lower-than- 



Red Bay is also associated with a region of 
average [O iii] surface brightness. 

Comparison of th e emission measure and electro n den- 
sity in this region (Garci'a-Di'az fc Henney 20071 shows 
that it is not a thin shell seen face-on, but has a thick- 
ness along the line of sight ~ 10^* cm, which is roughly 
twice its extent in the plane of the sky. We therefore 
suggest that the Red Bay may correspond to a stagnant 
zone of ionized gas (SI in Figure 4), which is trapped be- 
tween a concave ionization front, the stellar wind bubble, 
and the photocvaporation flows from molecular concen- 
trations such as the Bright Bar. 

3.2. The 3-D Structure of the Orion-S Region 

The several pieces of published and newly derived ev- 
idence require that we revise the accepted picture that 
Orion-S is simply a local rise in the MIF caused by dense 
underlying material in the OMC. Certainly, this is a re- 
gion of enhanced density as determined from H^'^CO^ 
observations. The recently published 21" reso lution 
H^'^CO+ maps that include the Huygens region (Ikeda 



et al.|2007 1 show the most massive cores to lie about 30" 
north of BN-KL and within Orion-S, but these observa- 



Garci'a-Di'az fc Henney] | |2007| , 

se [yinj observations, the Ked 



^ As originally drawn in Fig. 7 of < 
based on relatively low signal-to-noise 
Bay was broader, extending somewhat north of the Trapezium. 
However, this interpretation is not supported by the better-quality 
[O in] data. 
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tions cannot discriminate between cores that lie within 
a foreground optically thick feature and something that 
lies behind the MIF and within the main OMC filament. 
However, the presence of H2CO absorption in the radio 
continuum (discussed in § 2.3) clearly places a high den- 
sity, very optically thick object in front of an ionized vol- 
ume and this object could also be the source of the peak 
in rotationally excited CO emission that lies b etween the 
H2 CO ab sorption feature and the Trapezium (jWilson et 



al. 2001 1. The disparity in the extinction derived trom 
radio to optical and optical line ratios (discussed in § 2.4) 
further confirms that the optical features we see lie in the 
foreground of an ionized volume. The low ionization seen 
on the northeast side of Orion-S is most easily explained 
as an inclined face analogous to the escarpment causing 
the Bright Bar, but not as tilted. Otherwise it is impos- 
sible to explain the [Nil] surface brightness this close to 
the ionizing star. 

There are arguments that the foreground Orion-S 
Cloud hosts the embedded stars that are the sources of 
several optical features and outfiows. A dark rectangu- 
lar form lies almost exactly on the redshifted SiO outflow 
and the infrared source 139-357. It is not certain which 
(or if both) of these two features is associated, but an 
association of an optical feature with either places that 
feature in the foreground, rather than something within 
the OMC. There is a suggestive notch in the H2CO ab- 
sorption contour coincident with both the blueshifted 
outflows in SiO and CO, so that one or both of these 
must be in the foreground feature. The presence of high 
ionization, high tan gential velocity features symmetric 



with the D ark Arc ( jO'Dell fc Yusef-Zadeh|[2000| p'Dell 
et al.|2008 ) indicate that it is also part of the foreground 
structure. 

The H2CO absorption indicates very high optical 
depths in the near ultraviolet and it is unlikely that the 
obscuration by the foreground feature stops at the low- 
est contour of detection of H2CO absorption. This high 
extinction allows us to comment on the location of the 
sources that give rise to optical outflows. The sources 
of the large scale optical outflows have been isolated to 
a small region sho wn as the d ark dashed-line ellipse in 
Figure 1 CO'Dell 
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2008f. These outflows are all 

blueshifted! I'he east-oriented HH 529 shocks are high 
ionization and are seen very close to the source, indicat- 
ing that this flow emerges quickly from behind the local 
ionization front. The west-oriented HH 269 shocks and 
jet also emerge close to their sources and are seen visu- 
ally, arguing that they too originate in the foreground 
Orion- S Cloud. The emer gence further from the source 
region (O'Dell et al. [2008]) of the high velocity gas driv- 
ing the northwest-oriented HH 202 and the southeast- 
oriented HH 203 and HH 204 could either indicate that 
their sources are more deeply imbedded in the Orion-S 
Cloud or lie within the background OMC. The latter in- 
terpretation would imply that the lack of visibility near 
the source region is due to the high optical extinction of 
the Orion-S Cloud, but it seems more likely that all the 
sources lie at a similar distance along the line of sight. 

There is a question about the location of the strongest 
infrared and radio sources in this region: 137-408 (CS 3 
in earlier, low resolution studies) and 134-411 (FIR 4 in 
earlier, low resolution studies). The symmetry axis of the 
nearby HH 530 shocks points to either of these sources 



Observer 



-1 to 3 pc Veil 
~1.3 pc Veil-H+ 



Orion— S Bipolar Molecular 
Outflow Sources 
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Sources 



Fig. 5. — This sketch shows a southwest cross section passing 
through the Trapezium and Orion-S. The MIF is depicted as fiat 
behind Orion-S only because of a lack of information as to its true 
shape there. Although dimensions are shown, the uncertainty in 
the exact positions means that the figure is not drawn to scale. 

and the proximity of such an optical feature argues that 
its source lies in the foreground Orion-S cloud. 
The far-infrared emission from dust heated b y embed- 



ded strong infrared sources ( Houde et al.|2004 1 indicates 
(Figure 1) a peak in the vicinity of BN -KL and a weaker 
elongated peak in the Orion-S feature. The northerly ex- 
tension of this southern feature coincides with the source 
region for the optical outflows and the southerly exten- 
sion coincides with the 137-408 and 134-411 sources. Un- 
fortunately, this doesn't tell us if the southerly extension 
belongs to the foreground feature or lies within the OMC. 

The most likely geometry for this region is shown in 
Figure |5] In this model we have placed the Orion-S 
Cloud far enough into the foreground that the MIF ex- 
tends continuously behind it. The exact placement of the 
Orion-S Cloud in the foreground is uncertain, but it must 
be close to the distance of 9^ Ori C in front of the local 
MIF because if it were much further from the observer 
it could not be distinguished from the MIF. In the same 
way that the Bright Bar probably represents where the 
ionization front of the nebula encounters a long high- 
density condensation within the host OMC, it is likely 
that the Orion-S Cloud was a high density feature within 
the OMC, which has been isolated as the ionization front 
progressed into the host molecular cloud. This conclu- 
sion is compatible with the fact that the velocity of the 
H2CO extinction feature (Vh2CO=28.6±2) is essentially 
the same as the host molecular cloud (VbMC=25.9±1.7). 

3.3. Chemistry of the Orion-S Cloud 

Distinctly different heating processes are important in 
various regions of the Orion environment. Starlight is the 
most important heating/ionization process in the bright 
H+ region and its associated PDR. Wind-driven shocks 
are responsible for much of the geometry and emission 
of the BN object. Cosmic rays and heating by dissi- 
pative MHD waves may be important in deeper regions 
of OMCl. Given the assortment of possible dynamical 
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Fig. 6. — H2CO evaporation rate as a function of grain temper- 
ature. 

and microphysical processes involved, two approaches to 
numerical simulations of the ISM can be taken. Hydro- 
dynamics codes follow the motions of the gas but must 
compromise the microphysics to solve the problem on 
today's computers. A microphysics code (like Cloudy, 
which we use here) will not compromise on the micro- 
physics but must do something simple for the structure 
of the cloud. Atomic processes are usually much faster 
than dynamical timescales so the physics that determines 
the spectrum is often not affected by the dynamics and 
a static geometry can be assumed. 

In this section, we wish to understand several overall 
features of the Orion-S cloud via a simulation of the ion- 
ization, thermal, and chemical state of the cloud. The 
primary question we want to answer is "Can we explain 
the presence of H2CO in absorption?" . Furthermore, can 
we do so in the context of other molecules seen in emis- 
sion and the geometrical constraints imposed by the pro- 
jected size of Orion-S? There are estimates of the gas 
column density, particle density, temperature, and gas 
pressure (which is surprisingly high) . Can starlight from 
the Trapezium cluster account for these general proper- 
ties, or are other processes at play? Here we test whether 
starlight from the Trapezium cluster can account for sev- 
eral of the observed properties of Orion-S. 

3.3.1. Observational Constraints & Geometry 

The most important constraints to modeling the 
Orion-S cloud are the H2CO absorption measurements, 
the total column density, and the projected size of Orion- 
S. In cold, dark, molecular clouds, chemical models which 
include the effects of accretion of molecules onto grain 



surfaces ( Hasegawa et al.|1992||Hasegawa fc Herbst|1993| ) 
show that H2CO will reside almost exclusively on the 
surface of grains (henceforth "freeze out") at temper- 
atures at or below 10 K. A temperature of 10 K is 
also consistent with the ionized and molecular gas in 
Orion-S being in gas pressure equilibrium (see below). 
However, since H2CO is observed in the gas phase in 
Orion-S, the temperature must be higher, such that 
the rate of evaporation of H2CO from grain surfaces, 
feevap = 1-7 X 10^^ exp(— 1760/rd), must be large enough 
to produce a significant gas-phase abundance of H2CO. 
Figure [g] shows that this rate will be fast (> 10^^" s^^) 
for Td > 34 K. Therefore, H2CO observations imply 



that the dust temperature (which should be nearly in 
equilibrium with the gas temperature for deeper regions 
of the cloud) needs to be roughly 30 K or higher in or- 
der to explain the presence of H2CO. The total H col- 
umn density is given in McMuUin et al. (1993), who 
through 1.3 and 3.1 mm dust continuum emission found 
N(H) - 4 X 10^4 cm^2 f^j. region 12" N of the 3.1 mm 
continuum peak, corresponding to a region they referred 
to as the local quiescient cloud. The projected separation 
of the Trapezium stars from Orion-S is 2.2 x 10^^ cm, and 
the projected diameter of the cloud is ~ 2.8 x 10^^ cm. 
Assuming an approximately sp herical shape then give s a 
total cloud mass of (280 M©) ( [McMullin et al ||1993| . 

The close proximity of Orion-S to the Trapezium 
means there must be an H+ region on the side of Orion- 
S facing the Trapezium. Optical emission lines sug- 
gest that it has a density of ^ 600 cm~^ (Pogge et 



ar]|1992 Garcia-Diaz fc Henney||200"7 1 , corresponding to 
a gas pressure of F/k = 2n'l = 1.08 x 10^ cm~^ K, 
assuming a gas temperature of 8500 K. For the as- 
sumed thickness of Orion-S, and the above estimate of 
A^(H), the density deep in the molecular cloud of Orion- 
S is 3 X 10^ cm~^. If the environment were in gas- 
pressure equilibrium, then the corresponding molecular 
cloud pressure, 0.5nT = 2 x lO^T cm^'^ K would cor- 
respond to a temperature of 7 K. However, as noted 
above, this is inconsistent with the presence of gas-phase 
H2CO, indicating that the gas pressure in the molecular 
gas must be several times higher than that in the ionized 
gas. 

There are at least four different factors that may ac- 
count for this overpressure of the molecular gas: 

1. If the ionization front on the surface of the cloud is 

D-critical, then the transonic ionized flow from the 
front contributes an additional ram pressure term, 
which is roughly equal to the thermal pressure. On 
the neutral/molecular side, the net flow towards 
the front is highly subsonic, so the ram pressure 
term there is negligible. In a steady-state equilib- 
rium, this gives a molecular gas pressure twice that 
of the ionized gas. 

2. Absorption of the momentum of the stellar photons 
produces an outwardly directed effective gravity. In 
a steady-state conflguration, in which the molecu- 
lar gas is forced to remain at rest with respect to 
the star, this must be balanced by a pressure gra- 
dient between the optically thin and optically thick 
portions of the cloud. 

3. If the cloud is self-gravitating and hydrostatic, then 
the deeper gravitational potential in the center of 
the cloud will produce a higher pressure for the 
molecular gas there than for the ionized gas at the 
cloud surface. 

4. If the cloud is a transient phenomenon, formed 
from the dynamic interaction of converging flows, 
then there is no reason to expect it to be in pressure 
equilibrium at all. 

On the other hand, there are are two factors that would 
tend to work in the other direction: 
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Magneti c fields tend to be r elatively weak in ionized 
nebulae (Heiles et al. 1981 1, with magnetic pressure 
smaller than the gas pressure. In molecular clouds, 
on the other hand the reverse is frequently true, 
with the magnetic pressure typically exceeding the 
gas p ressure by a factor of 10 or more ( jCrutcherl 



1999) 



6. Turbulent internal motions in molecular clouds are 
frequently supersonic, with turbulent ram pressure 
approximately equal to the magnetic pressure. 

Therefore, the increase in gas pressure, when going from 
ionized to molecular gas, indicates that one or more of 
the first four factors must outweigh factors five and six. 

3.3.2. An Equilibrium Cloudy model of Orion-S 

We perform calculations with version 8 .0 of the spec- 
tral sy nthesis code Cloudy, described by [Ferland et al.| 
(19981. Exte nsions to PDR physics a nd chemistry are 



described by Abel et al. (2005 2008). The gas-phase 



composition and dust properties of the Orion environ 



ment are well known and are given in Abel et al. (2004a I. 

Our choice of ionizing continuum is characteristic 
of the dominant source of ionizing photons, the 07 
star 9^ Ori C. We assume a Kurucz stellar atmosphere 
with a temperature of 39,600 K, and Q{H) = 1.7 x 
10*^ photons s~^. We include X-ray emission from 
the stellar wind (see § 3.1), modeled as a 10® K 
bremsstrahlung component, although this does not af- 
fect the results presented below. 

For simplicity, we consider a time- independent, mag- 
netostatic equilbrium for the Orion-S cloud, as has been 
recently employed to accoun t for the structure an d ob- 
served magnetic field in M17 (Pellegrini et al.|2007[ ). The 
current version of Cloudy does not treat self-gravity, and 
we also ignore gas motions and time-dependent effects, 
so the model cannot directly address factors 1, 2, or 4, 
listed in the previous section. In the model, the total 
pressure increases as one proceeds into the cloud in such 
a way that the local pressure gradient exactly counter- 
balances the acceleration due to absorption of radiation. 
The total pressure is the sum of the thermal, magnetic, 
turbulent, and trapped resonance line radiation contri- 
butions. For simplicity we assume that the turbulent 
pressure is in equipartition with the magnetic field and 
that the magnetic field is tangled on small scales and well 
coupled to the gas, so that the magnetic stress tensor re- 
duces to an isotr opic pressure that is polytropic, with 



index 7„i = 4/3 (Henney et al. 2005). The magnetic 



field of the entire cloud is then set by the field at the 
illuminated face, increasing as magnetic 
field in Orion-S is not known but fi elds are importan t 
elsewhere in the Orion environment ( Abel et al.||2004a l. 
The initial magnetic field was taken as 3 /iC, which rises 
to 850 fiG deep inside the cloud, where the combined 
magnetic and turbulent pressures become roughly equal 
to the thermal pressure. We found that models with 
stronger magnetic fields were not capable of reproduc- 
ing the observed molecular hydrogen density deep i n the 



cloud of n ~ 3 X 10 cm since factors 5 and 6 of Si 3.3.1 



were too effective in counteracting factor 3. 

We assume a closed geometry, in which all radiation 
must eventually escape in the direction away from the 



central star cluster, which is appropriate if optically thick 
gas fully surrounds the cluster. Under this assump- 
tion, inwardly directed diffuse emission is compensated 
by symmetric emission from the other side of the star 
cluster. Light can undergo multiple scatterings as it dif- 
fuses outward though the cloud, so that the total rate of 
momentum transfer exceeds the single scattering limit, 
Lj./c, by a factor of a few. We use the iMcMullin et al. 



( 1993 1 deduced hydrogen column density to set the total 



thickness of the cloud. 

The entire structure and properties of the cloud follows 
from these assumptions. Cloudy computes the kinetic 
temperature, ionization, chemistry, optical depths, and 
emission properties of the gas at each point. The den- 
sity and radiation field striking the illuminated face are 
set by observations and the gas pressure at this point is 
computed. The calculation follows the starlight as it pen- 
etrates the cloud. The transmitted continuum sets the 
pressure and physical conditions. The magnetic field fol- 
lows from our flux- freezing assumption. The total pres- 
sure is set by the assumption of magnetostatic equilib- 
rium. The calculation stops when the observed column 
density is reached. 

Figure [7] shows the results of our calculations for the 
density, temperature, molecular abundances, and pres- 
sure as a function of depth into the cloud. The first crite- 
rion we demand our model reproduce is clearly satisfied: 
that the model yield a temperature high enough to allow 
H2CO to evaporate off of grain surfaces and remain in 
the gas phase. The dust temperature (Figure[7^) remains 
above 38 K thro ughout the model, which, as discussed 
above in § 3.3.1 is high enough to inhibit freeze-out of 



H2CO. This is achieved in our model purely through 
external heating by starlight from the Trapezium. 

Although we can explain the presence of gas-phase 
H2CO using our model, the chemistry of Orion-S is ex- 
tremely complex. At our molecular cloud temperature 
of 38 K, some molecules, such as SiO, CS, and SO2, will 
freeze out onto grains. Figure]?):; shows the predicted CO, 
CN, OH, and SiO abundance as a function of Ay- Our 
prediction of SiO, once the full effects of freeze out and 
molecular evaporation off grain surfaces is included, will 
likely be quite different. The current version of Cloudy 
only treats freeze-out and evaporation processes involv- 
ing CO, OH, and H2O, with the goal of making accu- 
rate predictions of gas-phase oxygen deep in molecular 
clouds. The freeze-out of other atoms/molecules, such 
as SiO, is not considered by our models, and neither are 
grain-surface reactions, other than those that form H2. 
So, while our current analysis is sufficient to deduce that 
H2CO remains in the gas phase, a full treatment of gas- 
grain and time-dependent chemical effects is necessary 
to fully understand all the rich chemical environment of 
Orion-S. ^^CO and ^"^CO lines are optically thick and 
not sensitive to gas properties such as pressure or col- 
umn density. C^^O and C^®0 should be optically thin 
and therefore potentially a good diagnostic of physical 
properties. However, as Orion-S and OMC-1 lie along 
the same sightline with similar kinematic properties, the 
optically thin CO emission from OMC-1 and Orion-S is 
difficult to disentangle. Therefore, we do not present CO 
comparisons because of these complexities. 

While our current analysis is sufficient to deduce H2CO 
remains in the gas phase, a full treatment of gas-grain 
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Fig. 7. — Physical structure for our constant pressure model of the ionized, atomic, and molecular gas in the Orion-S cloud. Plotted 
as a function of depth into the model are the (a) density, (b) gas and dust temperatures (c), molecular abundances, and (d) pressure 
constituents. 

and time-dependent chemical effects is necessary to fully 
understand the chemical environment of Orion-S. Such 
a treatment, while beyond the scope of this work and the 
current version of Cloudy (Cloudy currently only treats 
freeze-out of CO, OH, and H2O), will be the subject of 
a future investigation. 

One appealing aspect of our model is how we succeed 
in reproducing the projected dimensions of Orion-S un- 
der the assumption of magnetostatic equilibrium, with 
an effective gravity provided by the momentum of the 
absorbed starlight. The density (Figure [Tja) varies by 

4 dex over the physical extent of our model, yet our 
calculation reaches A^(H) ~ 4 x lO^'* cm~^ at a depth of 
2.4 X 10^^ cm, or within 20% of the observed projected 
thickness. Despite this success, there are arguably some 
deficiencies in our model assumptions, of which the most 
serious are the neglect of self-gravity and the rather weak 
magneti c field . Using the observed mass and size of the 
cloud (§ 3.3.1), one finds the escape speed at the cloud 
surface to be 7.3 km s~^, meaning that the molecular gas 



is gravitationally bound, while the ionized gas is gravi- 
tationally unbound. Thus, factor 3 of § |3.3.1 is likely to 
contribute significantly to the overpressure of the molec- 
ular gas, so that a self-gravitating model would allow 
the magnetic field in the cloud to be significantly higher, 
while still satisfying the observational constraints. How- 
ever, it is not clear that any equilibrium self-gravitating 
model is possible. Unless the turbulent and magnetic 
pressures are more than 100 times greater than given 
by our model, the cloud mass is higher than both the 
Jeans critical value and the magnetic critical value, and 
so the cloud may be dynamically collapsing on a free-fall 
timescale of a few times 10'* years. 

In summary, our equilibrium Cloudy model represents 
a necessary compromise, including as it does a detailed 
simulation of the relevant microphysics, but at the ex- 
pense of a simplified treatment of some of the relevant 
macrophysics. As such, although it is broadly successful 
in accounting for the observed properties of the Orion- 
S cloud using almost no free parameters, it should not 
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be considered the last word in modeling this region, es- 
pecially in the light of the uncertainties outlined in the 
previous paragraph. 

3.4. What is the Impact of a Putative High Relative 
Velocity of 9^ Ori C? 

In § 2.7 we saw that 9^ Ori C appears to be moving 
away from the OMC at about 13 km s~^. Since this star 
is only about 0.2 pc in front of the local MIF, 9^ Ori C 
would have been at the present location of the MIF only 
15,000 years ago. This is simply a reference number for 
the dynamic timescale as one would expect that 9^ Ori C 
would have continuously ionized the gas immediately sur- 
rounding it, thus producing a constantly changing struc- 
ture. It does appear that 9^ Ori C has only recently en- 
tered a region of low density gas and has just begun to 
photo- ionize a much larger volume of space (the Huygens 
region and the EON). This is a dramatically new way 
of viewing the Orion Nebula and means that the object 
need not be in quasi dynamic equilibrium. A short illu- 
mination age of the ONC stars is quite attractive because 
the well determined mass loss rates of the proplyds close 
to 9^ Ori C indicate that they should have been destroyed 
through photo-ablation withi n the last 1.5x10^-1.1x10^ 
years (Churchwell et al. 1987 | |Henney et al.|2002| [Smith 
erar|r3005[ |0'Dell & Henney '2008 i. The absence ot any 
evideiace for destruction indicates that their illumination 
ages by 9^ Ori C are no more than a few times 1.5x10"^ 
years. 

If one accepts the scenario that 9^ Ori C has emerged 
from a much higher density environment during the last 
10,000 years, then it is worth considering how this pic- 
ture would change the expected ionization, density, and 
dynamics of the present-day nebula. If the surrounding 
density was much higher in the past, then the Hii re- 
gion would have been much smaller. For example, in 
an ambient density of n — 10^ cm~'^, the Stromgren 
radius would be only 0.007 pc, or 3.2". As 9^ Ori C left 
the high-density environment, an R-type ionization front 
would propagate into the low-density gas at a speed of 
^n/n > 1 000 km s~^, where <I>h is the ionizing flux, de- 
fined in § |2.6.2| above. This speed is much larger than 
the stellar velocity, so that deviations from static pho- 
toionization equilibrium would be fleeting and difficult 
to detect. Assuming a density n — 100 cm~'^ for the low- 
density gas, the ionization front would reach its equilib- 
rium Stromgren value (~ 3 pc) in only 1000 years. The 
maximum recombination front velocity is also of order 
1000 km s~^ so that the gas would recombine effectively 
instantaneously in the high-density regions vacated by 
the star. The timescale for molecular hydrogen formation 
on dust grains is approxima tely 1000 (10^ cm^^/n) yr 
(Tielens & HoUenbach 19851, so that, unless the previ- 
ously ionized gas was of rather low density, there has 
been time for it to return to a fully molecular state. 

The density structure of the ionized gas is intimately 
tied to its dynamics, with typical internal velocities of 
order the ionized sound speed, which is very similar to 
the putative stellar velocity of 13 km s"'^. It seems likely, 
therefore, that the ionized density structure will not be 
greatly affecte d by the stellar motion. Numerical simula- 



tablished on a timescale equal to the sound crossing time. 
It is therefore at the largest scales in the nebula (> 1 pc) 
that one would expect to find the clearest evidence for 
non-equilibrium dynamics induced by the stellar motion. 
However, uncertainties in the exact ambient density dis- 
tribution and the effect of the stellar wind would make 
it difficult to draw clear conclusions. Another possible 
observational diagnostic of a moving ionized star would 
be the dense neutral wake, which numerical simulations 
predict shoul d form behind the star (Raga e t al.||1997 



Kraus et al. 20071. Again, observational detection o: 
such a wake would be very challenging, chiefly due to 
confusion brought on by the rich profusion of structures 
at all scales in the molecular cloud behind the nebula. 

A very young age for the present configuration of the 
Orion Nebula could also remove the dilemma presented 
by the fact that the H+ layer lying between 9^ Ori C 
and the Veil is moving towards the Veil with a closure 
velocity of about 18 km s~^ a nd there is no evid ence of 
an interaction with the Veil. Abel et alj ( 2006 ) calcu- 
late that the collision time with the Veil is about 40,000 
years. Without a detailed model of the expected very 
time-dependent structure, one cannot assess the impor- 
tance of this number except that it too indicates that the 
present configuration is quite young. 

In conclusion, it seems impossible to definitively con- 
firm or rule out the moving star hypothesis, based on ex- 
isting observational evidence. The strongest argument in 
favor remains the short mass- loss lifetimes derived for the 
proplyds in the inner nebula. The strongest arguments 
against are: (1) the required coincidence that 9^ Ori C 
currently resides so close to the center of the Trapez- 
ium cluster, despite not being dynamically bound to the 
cluster, and (2) the difficulty in explaining how 6^ Ori C 
acquired such a high velocity. The r nass of the cu rrent 
9^ Ori C multiple system is ~ 50 Mq ( |Kraus et alTp OOT), 
so assuming that it originally had a velocity typical of 
ONC members, then to achieve its current velocity it 
must have ejected a ~ 10 Mq companion at 65 km s"'^, 
or even faster for a less massive star. Note that these 
high stellar velocities mean that many-body interactions 
at the level of the star cluster are unlikely to play a role 
since the velocity dispersion of the cluster is only a few 
km s^^ (Jones & Walker 1988) and disruption of binaries 
by the cluster potential has been shown to be ineffective 
for binary separations < 0.5" (Reipurth et al. 2007). 

One further possibility that bears examination is that 
the velocity kick of the ionizing star was acquired much 
more recently. This would sidestep most of the argu- 
ments discussed above, both fo r and agains t the moving 



star hypothesis. For example, Tan (2004); Tan (20081 



has proposed that the ru naway H staFknown as the 



tions confirm ( Henney et al. 



2005a Arthur & Hoare 2006 1 



that an approximately steady-state champagne flow is es 



BeckHn-Neugebauer object ( Becklin & Neugebauer 1967 1 
was ejected f rom the 9^ Ori C multiple system 4500 years 
ago (but see Gomez et al. 2005[ '2008 for an alternative 
viewpoint). If the anomalous velocity of 9^ Ori C is due 
to this event, then it will only have moved about 0.05 pc 
since then, which is very small compared with the size 
of the nebula, so that observable consequences are un- 
likely. However, one serious problem with this proposal 
is that the measured radial velocity of BN is -1-13 km s~ ^ 
with respect to the ONC rest frame ( Scoville et al.|1983 l, 
which is the same absolute value as that recently deter- 
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mined for 6 Ori C ( Stahl et al. 2008 1, but in th e opposite 
direction . BN's mid-lK luminosity of 2600 Lq ( Gezari et 
al.|1998[ re-adjusted to a distance of 440 pc) is consistent 
with a dust-enshrouded main-sequence B2 to B3 star, im- 
plying a mass of 9 ± 1 Mq , which is more than 5 times 
less than the mass of 9^ Ori C. Therefore, as discussed in 
the previous paragraph, momentum conservation means 
that BN ought to have a much higher redshifted veloc- 
ity than is observed if it has truly been ejected from the 
9^ Ori C system in the recent past. 

3.5. Summary and Conclusions 

In this paper we have been able to establish several 
remarkable new properties of the Orion Nebula. 

1. There is an incomplete shell of high ionization gas 
shaped by the high velocity wind arising from 9^ Ori C 
and open in the direction of detected x-ray emission. 

2. The highest surface brightness portion of the Orion 
Nebula occurs to the southwest of the Trapezium stars, 
in the region called Orion-S, a center of star formation 
and stellar outflow. 

3. The Orion-S region is caused by an optically thick 
molecular cloud which lies within the ionized cavity of 
the Orion Nebula. This Orion-S Cloud is the host of a 



separate region of star formation and is sufficiently cold 
in the middle that it can produce H2CO lines in absorp- 
tion. 

4. It is possible that 9^ Ori C has a large radial ve- 
locity with respect to the ONC and other members of 
the Trapezium. The observations of the nebula do not 
contradict such a possibility and the presence of a large 
radial velocity would resolve the conundrum of the lack 
of destruction of the disks in proplyds close to 9^ Ori C 
and would explain why the inner layer expanding to- 
wards the observer has not yet reached the Veil feature. 
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